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LONG-TERM  GOALS 

One  long-term  goal  of  this  research  program  is  to  develop  an  operational  ionospheric  electron  density 
tomography  processor  with  the  following  characteristics: 

1 .  Solve  for  electron  density  on  a  global  grid  and  adapt  the  resolution  of  the  grid  based  on  the  data 
coverage 

2.  Be  capable  of  processing  data  over  a  long  time  period,  with  either  regular  or  varying  sample 
times  based  on  data  availability. 

3.  Accept  multiple  sources  of  data  from  both  ground  and  space  based  sources 

4.  The  output  electron  density  can  be  directly  ingested  into  ionospheric  data  assimilation  models. 

5.  Be  able  to  accept  as  input,  ionospheric  data  assimilative  model  inputs  to  initialize  the 
tomography  inversion. 

6.  Be  capable  of  running  in  near  real-time  in  a  completely  automated  manner 

The  last  year  of  research  has  resulted  in  an  initial  version  of  our  tomography  processor,  ionospheric 
data  assimilation  3D  (IDA3D),  which  has  partially  fulfilled  the  long-term  goals  stated  above. 

A  second  long-term  goal  of  this  research  has  been  to  validate  our  tomography  inversion  algorithms  in 
the  field  with  real  data  and  independent  data  from  incoherent  scatter  radars  (ISR)  and  ionosondes.  Our 
work  in  the  last  year,  as  well  as  research  conducted  by  other  tomography  investigators  such  as  the 
research  group  from  the  University  of  Wales  and  from  Northwest  Research  Associates  has  essentially 
demonstrated  the  validity  of  tomographic  mapping  of  the  ionosphere  in  mid-latitudes,  low-latitudes 
and  high-latitudes.  We  feel  now,  that  this  long-term  goal  has  been  met. 

The  third  long-term  goal  of  this  research  is  to  apply  our  tomography  processor  results  to  scientific 
analysis,  investigation  and  interpretation  of  the  dynamic  ionosphere.  During  the  last  two  years  we  have 
conducted  experiments  in  the  Caribbean,  Greenland  and  Alaska.  We  have  contributed  several 
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presentations  and  scientific  papers  on  the  interpretation  of  the  dynamic  ionosphere  from  tomographic 
results,  and  we  anticipate  continuing  that  work  over  the  next  three  years. 

OBJECTIVES 

The  first  technical  objective  of  this  research  was  to  develop  and  test  a  regional  three-dimensional 
tomography  algorithm  that  would  be  suitable  for  data  assimilative  models.  The  requirements  on  the 
algorithm  were  that  it  be  able  to  accept  beacon  satellite  TEC  data  as  well  as  ground  GPS  TEC  data  and 
ionosonde  data,  and  that  it  use  as  little  ad-hoc  a  priori  information  as  possible. 

Because  we  required  the  algorithm  to  be  suitable  for  data  assimilative  models,  it  needed  to  be  able  to 
make  an  estimate  of  the  accuracy  of  the  electron  densities  obtained  from  the  inversion.  In  addition, 
since  the  algorithm  would  apply  over  a  three-dimensional  region,  it  was  necessary  that  where  data  did 
not  provide  information,  that  the  output  of  the  algorithm  be  self-consistent  with  the  data  assimilative 
model  inputs. 

Einally,  in  order  to  meet  our  long-term  goals,  the  algorithm  had  to  be  designed  in  such  a  manner  that  it 
could  be  extended  to  a  global  grid  and  to  other  data  sources. 

The  second  technical  objective  was  to  field  tomography  receivers,  collect  data,  process  it  with  the 
tomography  algorithm  and  perform  scientific  analysis  on  the  results. 

APPROACH 

The  approach  adopted  here,  was  to  borrow  from  the  research  and  development  already  conducted  over 
many  years  by  the  atmospheric  weather  community.  In  particular,  the  concept  of  an  “objective 
analysis”  on  a  set  of  data,  resulting  in  an  “analyzed  state  variable”  was  adopted.  In  essence,  adopting 
this  overall  concept  allowed  us  to  generalize  beyond  the  concept  of  “tomographic  inversion 
algorithms”  to  the  concept  of  providing  the  very  best  analysis  of  a  model  state  variable,  given  an  initial 
model  prediction  and  a  set  of  data.  In  addition,  adopting  the  objective  analysis  framework,  and 
specifically  the  framework  of  3DVAR  articulated  by  NRE  West  ([1],[2])  allowed  us  to  make  us  of  an 
extensive  set  of  existing  literature  regarding  the  fusing  of  different  data  types,  the  effects  of  data  types 
on  the  analyzed  state  variables  and  given  the  analysis  the  re-initialization  of  all  of  the  model  state 
variables  and  the  next  cycle  of  model  predictions. 

The  algorithm  we  have  based  IDA3D  on  is  the  3DVAR  algorithm  of  NRE  West.  These  algorithms  are 
statistical  minimization  algorithms  (optimal  minimization  if  the  covariance  matrices  are  known 
exactly),  and  one  of  the  benefits  of  such  algorithms  is  they  provide  a  formal  measured  of  the  analyzed 
error,  given  the  initial  model  prediction  error. 

In  terms  of  meeting  the  scientific  analysis  objective,  coherent  ionospheric  Doppler  receivers  (CIDRs) 
developed  here  at  ARE:UT  have  been  deployed  to  the  Caribbean  for  two  combined  ionospheric 
campaigns  (CIC)  and  to  Greenland  and  Alaska  for  extended,  multiple-year  long-term  monitoring 
campaigns. 
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WORK  COMPLETED 


The  IDA3D  algorithm  has  been  developed  in  the  IDL  language,  tested  and  used  on  real  data  from 
deployed  arrays.  Following  is  a  list  of  the  current  tested  capabilities  of  IDA3D: 

1 .  A  3D  regional  uniform  or  regular  grid.  The  user  can  input  separate  arrays  of  latitudes, 
longitudes  and  altitudes,  or  can  specify  the  start  and  end  points  and  spacing  for  each  dimension. 

2.  The  algorithm  can  run  on  a  single  set  of  data  for  a  given  time,  or  it  can  be  run  sequentially  over 
many  days  worth  of  data. 

3.  The  algorithm  can  use  previous  results  to  initialize  the  next  inversion  using  a  Kalman  filter 
approach,  or  it  can  process  each  inversion  separately  using  a  climatalogical  model  such  as 
IRKS  [3]  as  the  initialization. 

4.  The  algorithm  accepts  beacon  TEC  data,  GPS  ground  TEC  data,  ionosonde  and  ISR  data  as 
inputs. 

5.  The  algorithm  uses  a  Gaussian  correlation  model  for  spatial  correlations  of  the  background 
model. 

6.  The  algorithm  outputs  the  analyzed  electron  density  on  the  specified  grid  as  well  as  the 
analyzed  variance  on  the  electron  density. 

7.  The  algorithm  has  been  run  in  near  real-time  on  test  data  downloaded  from  the  Greenland  array. 

Processing,  inversion,  scientific  analysis  and  interpretation  has  been  completed  on  data  sets  from  the 
Caribbean,  Greenland  and  Alaska.  In  addition,  older  data  sets  from  mid-latitudes  investigating 
ionospheric  storms  have  been  applied  to  IDA3D  and  new  analysis  performed.  The  results  of  these 
investigations  have  been  presented  at  several  conferences  and  three  archival  papers  have  been 
published. 
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Figure  1:  Four  images  of  ionospheric  structure  in  the  Caribbean.  The  upper  left  shows  a  2D  image 
(latitude  versus  altitude)  of  electron  density  from  the  tomographic  processor.  The  upper  right  shows 
an  image  of  vertical  TEC  (longitude  versus  latitude)  for  the  same  data  set.  The  lower  left  shows  a 
plot  of  data  from  two  CIT  receivers  that  went  into  the  inversion,  while  the  lower  right  show  the  GPS 

slant  TEC  data  that  went  into  the  inversion. 


RESULTS 

1.  CICl:  For  this  experimental  campaign,  a  suite  of  instruments,  by  many  researchers  were 
deployed  in  and  around  Puerto  Rico.  We  contributed  tomography  receivers  in  several  locations, 
and  also  analyzed  a  large  set  of  GPS  stations  from  throughout  the  region.  The  primary  result  of 
interest  from  the  campaign  was  the  imaging  by  the  all  sky  camera  of  the  so-called  mid-latitude 
turbulent  depletions  [4].  Figure  1  above,  shows  results  from  our  tomographic  and  GPS  analysis. 

2.  CIC2:  Was  very  similar  to  CICl  in  design  and  analysis.  Our  primary  result  was  to  compare  our 
IDA3D  results;  using  both  beacon  TEC  and  GPS  TEC  with  results  from  the  Utah  State 
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University  (USU)  ionospheric  predictive  model  as  well  as  from  the  JPL  GPS  based  TEC  maps. 
Figures  2  and  show  comparisons  between  IDA3D  and  the  JPL  TEC  maps. 

3.  Initial  Greenland  Analysis:  Figure  4  shows  an  inversion  from  the  Greenland  array  with  an 
Auroral-E  feature.  The  feature  has  been  confirmed  by  an  ionosonde  at  Narsarsauq.  This  is 
important,  since  if  tomography  can  reliably  and  accurately  image  E-region  electron  densities,  it 
can  contribute  important  spatial  and  temporal  information  to  the  Joule  heating  at  high  latitudes. 
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Figure  2:  2D  image  (longitude  versus  latitude)  of  VTEC  from  IDA3D  processor 
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Figure  3:  Same  region  as  Figures  2  and  3.  VTEC  from  JPL  GPS  TEC  maps 
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Figure  4:  Greenland  Reconstruction  along  scanning  meridian  of  Sondestrom  ISR  Note  the  E-region 
structure  near  -600  km.  This  structure  was  confirmed  by  an  ionosonde  at  Narsarsuaq 


IMPACT/APPLICATIONS 

This  research  impacts  both  space  weather  modeling,  estimation  and  prediction  and  the  propagation  of 

RF  signals  through  the  upper  atmosphere.  Relevant  applications  include  RF  communication  and 

navigation  systems,  RF  geo-location  systems,  and  HF  or  VHF  radar  systems  that  can  be  effected  by  the 

ionosphere. 

TRANSITIONS 

None 

RELATED  PROJECTS 

•  NSF  Greenland:  This  NSF  project  focuses  on  using  the  Greenland  tomography  data  with  IDA3D  to 
investigate  E-region  electron  densities  and  the  related  Joule  heating  effects. 

•  Alaska  Tomography:  This  project  in  collaboration  with  the  Geophysical  Institute  at  the  University 
of  Alaska  Fairbanks  and  Northwest  Research  Associates  is  designed  for  long  term  monitoring  of 
the  Alaska  ionosphere  in  support  of  HAARP  activities  and  to  do  scientific  studies  of  neutral  wind 
coupling  to  the  ionosphere  in  the  Fairbanks  region. 

•  Millstone  Hill  Tomography  Array:  This  project  is  sponsored  by  AFOSR  through  the  DURIP 
program  and  will  result  in  the  deployment  of  5  new  CIDR  receivers  along  an  array  extending  from 
Block  Island  Mass,  north  into  Canada,  and  along  a  meridian  coincident  with  the  Millstone  Hill  ISR. 
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In  collaboration  with  Millstone  Hill  we  expect  to  develop  as  part  of  the  Millstone  Hill  standard 
scheduling  software,  a  real-time  scheduler  for  IDA3D. 

•  NIWA  Data  Assimilation:  This  is  currently  a  small  project,  that  we  hope  will  grow,  designed  to 
couple  IDA3D,  ICS  with  a  fairly  simple  data  assimilative  model.  The  result  will  support  specific 
RF  applications  of  interest  to  NIWA,  and  the  current  plan  is  to  use  the  SAMI2  [7]  publicly 
available  predictive  model  for  assimilation. 
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